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Abstract. We provide the full set of renormalization group functions for the renormalization 
of QCD in the minimal MOM scheme to four loops for the colour group SU(N C ). 



1 Introduction. 

The main renormalization scheme used in quantum field theory is the modified minimal sub- 
traction, MS, scheme introduced in [H[2]. It has many elegant features which can be exploited 
to determine the renormalization group functions to a very high loop order. One of these is 
that of calculability. Briefly, only the divergences with respect to the regulator are removed, 
together with a specific finite part, ln(47re~ 7 ) where 7 is the Euler-Mascheroni constant [2j. Or- 
dinarily for conventional perturbation theory calculations one uses dimensional regularization in 
d = 4 — 2e dimensions and e plays the role of the regulator. Given this one need only consider 
the underlying massless quantum field theory safe in the knowledge that in this renormalization 
scheme the divergences will be mass independent. Thus, since massless Feynman graphs are 
significantly easier to compute than massive ones, then one can extract the ultraviolet diver- 
gences to high loop order. Moreover, using the scheme in gauge theories with massless fields 
gauge symmetry is preserved, [I]. Thus MS has been established as the favoured scheme for 
many years. However, for certain problems it is not necessarily the best choice. For instance, 
in lattice gauge theory computations it is not practical to implement since, for example, it is 
expensive for calculating Green's functions involving derivatives. Instead physical schemes such 
as the modified regularization invariant (RI') scheme have been introduced, OH]. These are 
asymmetric in that the definition is related to the choice of momentum configuration of 3-point 
functions. Related types of physical schemes, but not motivated by lattice considerations, are 
the momentum subtraction schemes of [5] denoted by MOM. These differ from RI' schemes in 
that the 3-point function momentum configuration is completely symmetric. Hence they do not 
suffer from infrared issues as the configuration does not have exceptional momenta. Both RI' 
and MOM schemes differ from MS in that finite pieces are absorbed into the renormalization 
constants which therefore depend on external momentum scales. As a corollary they are more 
difficult to calculate in analytically to high loop order. Irrespective of which scheme one chooses 
to use for an analysis, through the structure of the renormalization group equation it is possible 
to relate results. Thus within perturbation theory one can compute the conversion functions 
which allow one to map, for example, the coupling constant defined in one scheme to that in 
another. The other parameters and renormalization group functions can equally be related by 
the same formalism. 

A more recent development has been the introduction of another variant within the RI' and 
MOM family of physical renormalization schemes, [6j. It is called the minimal MOM scheme 
and is motivated by a property of the ghost-gluon vertex of QCD in the Landau gauge. This 
property is the non-renormalization of the vertex, [7]. However, the scheme is an extension of 
the concept beyond this specific gauge in a way which preserves a definition of the coupling 
constant in terms of the ghost and gluon form factors, [6]. This effective running coupling 
constant has been the subject of intense interest in recent years due to interesting features at 
medium and low energies which were noted earlier in [8j. For instance, it is believed that there 
are dimension two deviations from the expected running when compared to pure perturbation 
theory. The most recent work, [9], appears to reaffirm this property. With the minimal MOM 
scheme the effective coupling constant is not only simple to define but from a practical point of 
view does not require full knowledge of the ghost-gluon vertex function as would be necessary 
in other schemes [6]. Indeed in spite of being a non-exceptional momentum configuration it is 
numerically harder to extract a clean signal from the lattice for a fully symmetric vertex such 
as in the MOM context. In developing the minimal MOM scheme, [6], the four loop QCD 
/3-function was determined for SU(N C ) not only in the Landau gauge but also for a particular 
formulation of a linear covariant gauge. The full set of renormalization group functions as far 
as they could be calculated were not given. Therefore, it is the purpose of this note to provide 



the wave function renormalization group functions to as far a loop order as is possible. For 
an arbitrary colour group this is to three loops for the wave functions and four loops for the 
j3- function and quark mass anomalous dimension. Though we will also provide the former to 
four loops for SU(N C ). We will do this in two ways. The first is the direct evaluation of all 
the three loop Green's functions where the minimal MOM renormalization scheme definition 
is implemented directly. The second is by construction of the associated conversion functions 
and use of the renormalization group equations. This will serve as a check on our computations 
and allow us to deduce four loop information. One of the reasons for the direct renormalization 
is that it provides a non-trivial independent check on the results of [6j. There the /3-function 
was adduced from known finite parts of Green's functions given in |10j . As a separate exercise 
we choose to work in a minor variation of the original minimal MOM scheme and that is to 
renormalize the gauge parameter, a, in the full ethos of the minimal MOM scheme. In [6] 
the renormalization of the gauge parameter is completely equivalent to that of the MS scheme. 
Though our results will equate for the Landau gauge. A similar issue for a arises in the RI' case 
[llj . Finally, given the interest in the behaviour of the effective coupling constant in the Landau 
gauge and its power law deviation, we will provide the minimal MOM anomalous dimension 
of the operator thought to be associated with the dimension two correction which is ^A^A 0,11 
where A a is the gluon. This is in order to allow one to perform a complete renormalization 
group running analysis in the minimal MOM scheme for such infrared problems. 

The paper is organized as follows. We recall the definition and properties of the minimal 
MOM scheme in section 2 before recording our results in the subsequent section. These include 
the renormalization group functions and the conversion functions for an arbitrary colour group. 
For MS renormalization group functions which are only known at four loops for SU(N C ), we 
provide the corresponding minimal MOM scheme results in section 4 together with those for the 
dimension two operator anomalous dimension. A conclusion is provided in section 5. 



2 Formalism. 

We begin by recalling the definition of the minimal MOM scheme, [6j. First, if we denote bare 
quantities in the QCD Lagrangian by the subscript , then in our notation the renormalization 
constants, Zj, are given by 

A*» = y/Z^A^ , c a Q = y%c a , Vo 7 = v^ 47 

a = Z~ l Z A a , m Q = mZ m , g Q = n € Z g g (2.1) 

where A a is the gluon, c a is the Faddeev-Popov ghost and ip l is the quark. The indices have 
the ranges 1 < a < Na, 1 < i < Np and 1 < / < Nf where Np and Na are the respective 
dimensions of the fundamental and adjoint representations of the colour group and Nf is the 
number of massive quarks each of the same mass m. The coupling constant is g and a is the 
gauge parameter of the linear covariant gauge. The Landau gauge corresponds to a = 0. We 
use the above definition of the renormalization of a to be consistent with |llj . Throughout 
we use dimensional regularization in d = 4 — 2e spacetime dimensions and the mass scale \x is 
introduced to ensure the coupling constant is dimensionless in d-dimensions. With these formal 
definitions of the renormalization constants they are then determined explicitly by specifying 
a scheme to absorb the infinities in the various 2 and 3-point functions of the theory. For 
instance, MS corresponds to removing only the poles in e together with a certain finite part at 
some subtraction point. 

For momentum subtraction schemes, denoted generally by MOM, the scheme is defined such 
that at the subtraction point the poles in e together with all the finite part are absorbed into the 



renormalization constant, [5j. For the QCD Lagrangian this produces several different schemes 
since there are several vertices which one can use to define the coupling constant renormalization. 
Choosing one, say, means that the remaining vertex functions are finite and consistent with the 
Slavnov- Taylor identities. The variation on this approach introduced in [6] is that the 2-point 
functions are renormalized using the MOM criterion of [5] but the 3-point vertices are treated 
differently. Specifically, to ease comparison with lattice analyses the completely symmetric 
subtraction point of [5j is not used. Instead the asymmetric point is used where the external 
momentum of an external leg is nullified. Moreover, partly motivated by the non-renormalization 
of the ghost-gluon vertex in the Landau gauge, [7], the coupling constant renormalization is 
defined by ensuring that this vertex renormalization constant is the same as the MS one. One 
benefit of this, [6], is that to define the scheme one only needs to know the vertex structure in 
the MS scheme which reduces work for non-perturbative applications. In our notation, (12.1ft . 
this corresponds to, [6j, 



Z MS Z MS Z MS = ^mMOM AmMOM^mMOM ^ 2) 

where mMOM denotes the minimal MOM scheme. Though, in this formal definition it is 
important to appreciate that the variables g and a on either side of the equation are in different 
schemes. We note that throughout our convention is that when a scheme is specified as a label 
on a quantity then it is a function of the parameters g and a in that scheme. With (|2.2p then all 
the renormalization constants of massless QCD are defined for the minimal MOM scheme, [6]. 
As noted earlier in [6] the gauge parameter renormalization was treated as an MS one rather 
than define it as the full MOM renormalization as used in [5]. Therefore, we will follow the 
approach of [5] here and have a minimal MOM a. From the practical point of view our results 
in the Landau gauge will be the same and differ only in the a dependent part. 



The procedure we have used is to apply the Mincer algorithm, [12] . to the massless QCD 
Lagrangian and compute all the 2-point functions as well as the ghost-gluon vertex at the 
asymmetric point. The quark mass anomalous dimension will be discussed later. This algorithm 
evaluates massless three loop 2-point functions to the finite part in dimensional regularization. 
It has been encoded, [13J, in the symbolic manipulation language Form, [14] , which is our main 
computational tool. The Feynman diagrams are generated by Qgraf, [15j, and the output 
converted into Form input notation. As all graphs are evaluated to the finite parts then we 
can extract the explicit renormalization constants in the minimal MOM. We do this first by 
renormalizing the 2-point functions before defining the coupling constant renormalization via 
(12.2ft . Then one proceeds to the next loop order. In using the definition of the coupling constant 
renormalization we have to relate the parameters between the schemes. For the gauge parameter 

this is given by 

Z MS 

"mMOM^) = zm MOM a MS^ ( 2 - 3 ) 

where we used the fact that we are in a linear covariant gauge which implies Z a = 1. To ensure 
a finite expression in e emerges the parameters within z™ have to be converted to their 

MS partners. This is achieved order by order in perturbation theory. We have determined these 
to three loops and, with a = <? 2 /(167r 2 ), found 
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a mMOM = fl + 
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a 
36 

2^2 , nooc:„2^»2 , oc/- A s(o\ _ /o2 , joi„/i2 



9a z C A + 18aC A + 169Ca - 80N f T F 

405a 3 C A - 486C(3)a 2 Cl + 2835a 2 C A + 3564C(3)aCl + 2421aC| 
- lU0aC A N f T F - 6318C(3)Cl + 76063^1 - 10368((3)C A N f T F 



50656C A N f T F + 20736C(3)C , FA r /T F - 23760C F N f T F + 6400NfT F 



a 3 
1296 



and 
mMOM 



+ 



-10692C(3)a 4 Cl + 8505C(5)a 4 C;| + 41067a 4 C;| - 59292C(3)a 3 Ci 

- 4860C(5)a 3 Ci + 293301a 3 Ci - 138024C(3)a 2 Ci - 85050C(5)a 2 C^ 
+ 1315035a 2 Ci + 46656( (3) a 2 C\N f T F - 322056a 2 C 2 A N f T F 

+ 3355020C(3)aCi - 860220C(5)aCi + 1277496aCi 

- 2115072C(3)aC|Ar / T F - 581760aClN f T F + 373248C (3) aC A C F N f T F 

- 427680aC A C F N f T F + 33l77QQ{3)aC A NfT F + 32256aC A N^T F 

- 6552900C(3)Ci - 1896615C(5)Ci + 42074947C^ - 4499712C(3)C^ N f T F 
+ 2488320C(5)C|iV f T F - 38975424(7^ T F + 15303168( (3) C A C F N f T F 

+ 3732480C(5)C^ C F N f T F - 23755968C ' A C F N f T F + 3068928C(3)C j4 iV / 2 T| 
+ 9209280C A NfT F + 4603392C F N f T F - 7464960C( 5) C F N f T f T F 



+ 1482624CpN f T F - 
+ 8065152C F iV f 2 T| 



6469632C(3)C F iY f 2 r| 



512000iY/T| 



46656 



+ 0(a 5 



(2.4) 



a + 



-9a 2 CU - 18aC A - 97C A + 80N f T F 



a. a 
36 



l8a A C\ - l8c?C 2 A + \9Qa 2 C\ - 320a 2 C A N f T F - 576C(3)aCi + 463aCl 
- 320aC A N f T F + 864C(3)Cl - 7143C 1 ! + 230AC{3)C A N f T F + A2A8C A N f T F 



aa 



18792a 4 C\ 
75951a 3 C\ 



4608((3)C F N f T F + 5280C F N f T F ] — 

288 

A86a e C A + 1944a 5 <7^ + 4212C(3)a 4 C^ - 5670C(5)a 4 Ci 

+ l2960a A C A N f T F + 48276C(3)a 3 Ci - 6480C(5)a 3 C^ 

- 8640a 3 C , i7V f T F - 52164C(3)a 2 (7i + 2916C(4)a 2 Cl + 124740C(5)a 2 C^ 
+ 92505a 2 Ci - 129600<C(3)a 2 C 2 A N f T F - 147288a 2 C 2 A N f T F 

+ 248832C(3)a 2 C A C F N f T F - 285120a 2 C A C F N f T F - 384Q^a 2 C A NfT F 

- 1303668C(3)aCi + 11664C(4)aCl + 447120C(5)aCl + 354807aCl 
+ 698U2C (3) aClN f T F - 3\233QaC 2 A N f T F + 248832C(3)aC j4 (7 F A? f T F 

- 285120aC A C F N f T F - 221184C(3)aC A iV / 2 Tj. + 55296aC A NJ : T F 
+ 2007504C(3)Ci + 8748C(4)Ci + 1138050C(5)C^ - 10221367(7^ 

+ 1505088((3)C A N f T F - 279936C(4)C^A? / T F - 1658880£(5)C A N f T F 
+ 9236488C|iV / :7> - 5156352C(3)Ca C F N f T F + 373248C(4)CaC> iV/Tr 

- 2488320C (5) C A C F N f T F + 929366AC A C F N f T F - 884736({3)C A Nj?T F 

- 13A3872C A NfT F - 3068928((3)C F N f T F + 4976640C(5)C| iV/Tp 



r2rr-i2 



2842368C ' F Nj?T F 



- 988416C F N f T F + 2101248C(3)C F iV/TJ, 
+ 0(a A ) 
where £(z) is the Riemann zeta function. The group Casimirs are defined by 



aa 
31104 



Tr TT 



T F 5' 



ah 



rpdrpQ 



C F I 



racd rbcd 



c A s 



all 



(2.5) 



(2.6) 



where T a are the generators of the colour group whose structure functions are f abc . In (|2.4p 
and (|2.5p the variables on the right hand side are in the MS scheme. For the Landau gauge it 



is easy to see that then the parameters coincide. We have checked that (J2.4J) agrees with the 
alternative definition of the mapping given in [6] based on the actual finite parts of the gluon 
and ghost 2-point functions after their MS renormalization. 

While we will perform a direct evaluation of the renormalization constants in the minimal 
MOM, there are several checks which will be carried out. One is to exploit properties of the 
renormalization group equation which allows one to map the anomalous dimensions deduced 
in each scheme via conversion functions which are denoted by Ci(a,a) where i will be a label 
corresponding to a field or a parameter. First, we will perform the explicit renormalization in 
the minimal MOM and deduce the anomalous dimensions directly. Then we will compute the 
conversion functions and from these construct the anomalous dimensions indirectly. Thus if we 
define the conversion functions by 

vmMOM ymMOM 

C mMOM (a,a) = ^=- , Cr M0M (a,a) = =*=- (2.7) 



where (ft € {A,c, ip}, then the minimal MOM renormalization group functions are given by 
/3 m (%MOM) a mMOM) = 



oMS[„ ^"mMOM 



'MS' q„ 



MS/ \ " m MOM 



+ a MS^« ( a MS' a MS 



da 



'^ a MS J MS^mMOM 



and 

mMOM 



(«mMOM) a mMOM) = I4, ( a MS' a Ms) 



,Mti / \ d 



+e MS (m)Tt=^r 0M (< 



■ 



f) n LLL ^* v°ms' a 'Ms; 

aa MS 

i „ ^MS / \ g , r mMOM/ \ 

^ U MS'« ^MS'^MSyp^, ^<t> \ "MS' "MS J 

u MS 



MS^mMOM 
(2.9) 



Here MS —> mMOM means that after computing the right hand side the expression will be a 
function of MS variables and these must therefore be converted to minimal MOM ones. The 
relations are given by inverting (J2.4J) and (|2.5p . One benefit of this formalism is that it can be 
exploited to produce the four loop anomalous dimensions and /3-function. The reason for this 
is that the three loop conversion functions give a four loop contribution to the minimal MOM 
anomalous dimensions and /3-function and as the MS versions of these are known, [16^ \TT\ \TE[ 
[El [201 EH E21 [231 [Ml [2S1 [26] , then the left hand side can be deduced at four loops. 

3 Results. 

We now formally record our results. The one loop expressions will be the same as the MS ones 
since that term is scheme independent. This includes the /3-function as we are using a mass 
dependent renormalization scheme and only in mass independent schemes is the two loop term 
scheme independent. Moreover, the /3-function will be gauge dependent for the same reason. 
Therefore, we havqj 

/3 mMOM (a,a) = - [UC A -4N f T F ]^ 



*A data file is attached which gives an electronic version of all our expressions. 



+ 



- 3a 3 C 2 A + Wa 2 Cl - 8a 2 C A N f T F + YiaC\ - 8aC A N f T F - 136C A 

a 3 
+ 80C A N f T F + A8C F N f T F ] — 

-l65a 4 C A + 2Aa 4 C\N f T F + 108C(3)q 3 C*1 - 189a 3 C A 

- U4a 3 C A N f T F - 468C(3)a 2 Cl + 2175a 2 Cl + 144C(3)a 2 C A N f T F 

- 1656a 2 C A N f T F - 864a 2 C A C F N f T F - 1188C(3)aCl + 3291aC A 

- 1776aC A N f T F - 1152aC A C F N f T F + 5148C(3)Ci - 38620(7^ 
+ 6576((3)C A N f T F + 3214AC%N f T F - 16896C (3) C A C F N f T F 

+ 20512C A C F N f T F - 3072((3)C A NfT F - UlQC A NfT F 

i a 4 

- 576C F N f T F + 6144C(3)C F iV / 2 r| - 5888C F Nj*T F — - 

864C(3)a 5 Cl - 3780C(5)a 5 CJ - 11745a 5 CJ + 1728a 5 C A N f T F 

+ 32472C(3)a 4 Cl + 4U0((5)a 4 C A - 81549a 4 C\ 

- 1440C (3) a 4 C A N f T F - 5040C (5) a 4 C7| N f T F + 7200a 4 C A N f T F 
+ 777Qa 4 C 2 A C F N f T F + 47052C(3)a 3 Cl + 19800((5)a 3 C A 

- 81873a 3 CJ + 18432C(3)a 3 C A N f T F + 1440C(5)a 3 C^iV/TV 

- 67752a 3 C 3 A N f T F - 7776a 3 C A C F N f T F - 397368C(3)a 2 Cl 
+ 152280C(5)a 2 Cl + 1028898a 2 CJ - 36576C(3)a 2 C 3 A N f T F 

- 1098936a 2 C 3 A N f T F + Q393mC,{3)a 2 C 2 A C F N f T F 

- 790272a 2 C A C F N f T F + 73728C(3)a 2 (7 2 1 iV / 2 r| + 133632a 2 C 2 A NfT F 
+ 20736a 2 C A C F N f T F - 221l84((3)a 2 C A C F NJ i T F 

+ 21 1968a 2 C A C F NfT F - 2400708C(3)aCJ + 987660C(5)aCl 
+ 1719423aCl + 1655712C(3)aC|AT / r F - 254880C (5) aC 3 A N f T F 

- 1817880aC A N f T F + 798336( (3) aC A C F N f T F 

- l030752aC A C F N f T F - 6l7472((3)aC A NJ?T F + 391488aC| NfTp 
+ 3lW4aC A C F N f T F - 331776C(3)aCA C F NfT F 

+ 3\79h2aC A C F NfT F + 98304C (3) aC^iV/T 3 ; - 24576aCU iV/Tl 
+ 5509416C(3)Cl + 3090780C(5)CJ - 22106704CJ 

- 1217376C(3)Cl N f T F - 5178960((5)C A N f T F + 23501280C|iV / T F 

- 7050240C(3)C|C F 7V / T F - 6082560C (5) C\C F N f T F 
+ l7477280C A C F N f T F - 1654272C(3)C'iiV / 2 T| 

+ 1474560C(5)CliV / 2 T 2 ; - 5719680C|iV / 2 T| 

- 7907328C(3)C j4 C|iV / Tp + 12165l20((5)C A C F N f T F 

- W7104C A C F N f T F + 4755456C(3)C A C F NfT F 

+ 2211840C(5)Ca C F NfTp - 1086l056C A C F NfTp 

+ 344064C(3)C*AiV / 3 r| + 22937SC A NfTp - A76928CpN f T F 

+ 3538944C (3) CpN^Tp - 4423680C(5)C| NfTp + 267264Cf, JV/r| 

d {4) 

- 884736C(3)C ' F NfTp + 1327104C F iV / 3 r| - 2433024C(3)-^ 



7(4) 



vW 



+ 92160-^ + 5750784C(3)-g4JVj> - 589824 



7V A 



AU 



J (4) 



AV 



7(4) 



1769472C(3)-££./V f 2 + 811008- 



i (4) 



4. roiiuuo— — JV f 

A^a / i\fr i 



10368 



+ 0(a 6 



(3.1) 



As the four loop MS /3-function was computed for an arbitrary colour group, [23], the general 
colour group Casimirs appear. In our notation they are defined by 



d 



(4) 



jabcd jabcd j\ 

Up Up , Upj^ 



7(4) 



jabcd jabcd 



jaocajaoc 
pp - - up up , u, FA - - Up U A , u AA 



d 



(4) 



jabcd jabcd 



and the totally symmetric rank 4 colour tensor is defined by, |27j , 

mbcd . _^T\, (rpO,rp(brpCrpd)\ 



6 

where the group generators are in the R representation. 



(3.2) 



(3.3) 



For the anomalous dimensions only the three loop MS expressions are known for an arbitrary 
colour group, [21]. Thus to the same order the minimal MOM expressions are 



mMOM/ \ 
lA {a, a) 



[3aC A - 13C A + 8N f T F ] 



6 



+ 



+ 



-6a 3 C A + Ylc?C\ - l6a 2 C A N f T F + !7aC A - 16aC A N f T F - 170C A 



+ !36C A N f T F + 96C F N f T F ] 



24 



-165a 4 C A + 24a 4 C A A f T, F + 54C(3)a 3 C A - l26a 3 C A - U4a 3 C A N f T F 

- 576C(3)a 2 C A + 1761a 2 C A + U4((3)a 2 C A N f T F - 1512a 2 [ C 2 A N f T F 

- 8Ma 2 C A C F N f T F - 774C(3)aC A + 102aC A - 288C{3)aC A N f T F 

- mQaC A N f T F - 1152aC A C F N f T F + 3456C(3)C A - 23032C A 
+ 6288C(3)C A A f 7> + 21320C A N f T F - 16896((3)C A C F N f T F 
+ 196A8C A C F N f T F - 3072((3)C A N f 2 Tp - 24QQC A NfTp 



Irj^l 



576CpN f T F + 61UC(3)C F NfTp-5888C F N/Tp — + O^) 



288 



mMOM 

Ic 



(a, a) 



[a - 3] — — + 



3a 2 C A - 3aC A - 34C A + 8N f T F 



C A a 2 
16 



+ 



54C(3)a 3 Cl - 45a 3 C A - 36C(3)a 2 C A + 216a 2 C A - 48a 2 C A A f T F 

+ 42C(3)aC A + 109aC A + 96C(3)aC A A f T F - 152aC A N f T F 
+ 564C(3)C A - 5196C A + 96C(3)(7 A A f T F + 3608C A N f T F + 288C F N f T F 
C A a 3 



2rp2 



QWNfTp 



^ MOM (a,a) 



f 

aC F a + C F 
+ 



192 



+ 0{a") 



3a 2 C A + 6aC A + 25C A - 6C F - 8N f T F 



C F u 2 



.,2/o2 



18({3)a- i C A + 27a d C A - 2Aa 6 C A C F - 90C(3)c^C A + 123a z C A 

- 36a 2 C A C F + 48a 2 C A N f Tp - 618C(3)aC A + 395aC A + 72aC A C F 
+ 192((3)aC A N f T F - UaC A N f T F - 1470C(3)C A + 3843C A 
+ 576C(3)C A C F - 1260C A C F + 384C(3)C A A f T F - 18A0C A N f T F 

C F u 3 



+ 72C F - mC F N f Tp + \28NfTp 



48 



+ 0(u 4 



(3.4) 



We have checked explicitly that the gauge parameter satisfies 

,mMOM/ n mMOM/ n 

(a, a) = - 1a K") 



7,, 



which is a check on our calculation. 



(3.5) 



Having provided the anomalous dimensions we have checked that they are completely repro- 
duced using the conversion function approach. The explicit forms of these functions are 



C g (a,a) 



+ 



+ 



a 
72 

1242a 2 C\ 



1 + -9a z C A - l8aC A - 169C A + 80N f T F 

243a 4 C A - 648a 3 C A + 1944C(3)a 2 C A - 1242a 2 C A - 4320a 2 C A N f T F 

- 14256C(3)aCl + 8568aC A - 2880aC A N f T F + 25272C(3)C A - 218569C A 

+ 4U72((3)C A N f T F + 12l504C A N f T F - 82944C (3) C F N f T F + 95040C F N f T F 

i a 2 

-6400N?T F 

f F \ 10368 

-3645a 6 C A + 21870a 5 C A + 33048C(3)a 4 C A - 68040C(5)a 4 C A 

- 183951a 4 C A + 97200a 4 C 2 A N f T F + 754272C(3)a 3 C A + 38880C(5)a 3 C A 

- 1501740a 3 C A - 155520a 3 C 2 A N f T F + 206064C(3)a 2 C A + 680400C(5)a 2 C A 

- 710235a 2 C\ - W26432( (3) a 2 C A N f T F - 1887840a 2 C 2 A N f T F 
+ 2239488C (3) a 2 C A C F N f T F - 2566080a 2 C A C F N f T F - 172800a 2 C A NfT F 

- 20977056C(3)aCl + 6881760C(5)aC A + 3407958aC A 
+ 11259648C(3)aC A iV / T F - 4231 29§aC A N f T F + U92992( (3) aC A C F N f T F 

- 17W720aC A C F N f T F - 2654208((3)aC A Nj i T F + 7787b2aC A NfT F 
+ 39610296C(3)C A + 15172920C(5)C A - 206477857(7^ 
+ 21036672C(3)C|iV / T F - 19906560({5)C A N f T F + 170325744C|AT / T F 

- 80372736C(3)C A C F iV / Tp - 29859840C(5)C A C F iV / T J F 
+ U186246AC A C F N f T F - 14598144C (3) C^iV/T 2 , - 28289280C A iV / 2 Tj, 

- 36827136((3)C F N f T F + 59719680((5)C F N f T F - 11860992C ' F N f T F 



+ 31850496C(3)C F iV f 2 r| 



j2rj-\2 



3n-t3 



A1711616C F N/T F + 512000iV/r| 



746496 



+ 0(a 4 



(3.6) 



C A {a,a) 



+ 



+ 



a 
36 



1 + 9a 2 C A + 18aC A + 97C A - 80N f T F 

810a 3 C A + 2430a 2 C A + 5184C(3)aC A + 2817aC A - 2880aC A N f T F 
- 7776C(3)C A + 83105C A - 20736C(3)C A AT / T F - 69272C A N f T F 

n 

+ 41472C (3) C F N f T F - 47520C F N f T F + 12800AT/TJ 



2592 



-12636C(3)a 4 (7 A + 17010C(5)a 4 C A + 64638a 4 (7 A - 51516C(3)a 3 C A 

+ 19440C(5)a 3 C A + 322947a 3 C A + 203148C(3)a 2 C A - 8748C(4)a 2 (7 A 

- 374220C(5)a 2 Ci + 1094553a 2 C A + 15552C(3)a 2 C 2 A N f T F 

- 303912a 2 C 2 A N f T F + 4636764C(3)aC A - 34992C(4)aC A 

- 1341360C(5)aC A + 1457685aC A - 3670272C(3)aC A iV / T F 

- 890064aC A N f T F + 746496C(3)aC A C F iV / T F - 855360aC A C F N f T F 



+ 663552C (3) aC^iV/Tj + Mol2aC A NfT F - 7531056C(3)C A 

- 26244C(4)C A - 3414150C(5)C A + 44961125C A - 7293888C(3)C A N f T F 

+ 839808C (4) C\N f T F + 4976640C(5)C A N f T F - 49928712C A N f T F 

+ 235U624({3)C A C F N f T F - 1119744C(4)C ' A C F N f T F 

+ 74M9GOC(5)C A C F N f T F - 37099872C ' A C F N f T F + 5971968({3)C A Nj?T F 

+ 13873536CaA^ 2 T| + 9206784C(3)C| N f T F - 14929920C(5)C F N f T F 



r2rri2 



+ 2965248C F N f T F - 12939264((3)C F NfT F + 16130304C F N/T F 



1024000iV f 3 r| 



93312 



+ 0(a 4 



(3.7) 



C c (a, a) 



l + C A a 



and 



C^(a,a) 



+ 



+ 



-36C(3)a 2 C A + 72a 2 C A + 72((3)aC A - 2\aC A 

3/^2 i ono/ii „3/^-<2 



- 180C(3)Ca + 1943C A 

-11178C(3)a 3 C A - 4860C(5)a 3 C A + 29241a 3 C A - 56862C(3)a 2 C A 

+ 1458C(4)a 2 Ci + 34020C(5)a 2 Ci + 102789a 2 C A + 254826C(3)aC A 
+ 5832C(4)aC A - 63180C(5)aC A - 3510aC A - 67392({3)aC A N f T F 
+ 42984aC A N f T F - 728082C(3)C A + 4374C(4)C A - 63180C(5)C A 
+ 4329412C A - 100224C(3)C A A r /T F - 139968C(4)CUiV/ 1> - 2650192C A A r /T F 
+ 684288C (3) OiV/Tp + 186624C(4)C F N f T F - 1165104C F N f T F 

C A a z , „, 4 



+ 27648C(3)iV / 2 T 2 , + 330304^/^ 



31104 



+ 0{a A ) 



(3.8) 



1 
+ 

+ 



52aC A + 24C(3)C A 



- aC F a 
-9a 2 C A + 8a 2 C F + 24C(3)aC A 

C a 2 

- 82C A + 5C F + 28N f T F ] — ^— 

8 

1728C(3)a 3 C A - 11880a 3 C A - 5184C(3)a 3 C A C F + 12312a 3 C A C F 

+ 3456C(3)a 3 C 2 ? - 5184a 3 Cf, + 25272C(3)a 2 C A - 972C(4)a 2 C A 

- 6480C(5)a 2 Cl - 63747a 2 C A - 31104C(3)a 2 C A C F + 59616a 2 C A C F 



+ 181440C(3)aCl - 1944aC(4)C A - 12960C(5)aC A - 358191aC A 
+ 57024C(3)aC A C F - 103680C (5) aC A C F + 85536aC A C F 

- 41472C(3)aC A iV / r F + 124056aC A N f T F - U016aC F - 28512aC F N f T F 
+ 678024C(3)C A + 22356C(4)C A - 213840C(5)C A - 1274056C A 

- 228096C(3)C A C F - 31104C(4)C A C F + 103680C(5)C A C F + 215352C A C F 

- 89856C(3)C A iV / r F + 760768C A N f T F + 31536C| - 829M( (3) C F N f T F 

C F a 3 



+ 68256C F N f T F 



100480 NfT F 



5184 



+ 0(a 4 



(3.9) 



We use the convention that the variables on the right hand side are in the MS scheme. 



While [6] provided the renormalization group functions for massless QCD it is possible to 
deduce the quark mass anomalous dimension to four loops for an arbitrary colour group. This 
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requires the conversion function for the quark mass renormalization and the four loop MS 
anomalous dimension. The latter has been provided in [23] and [25]. To deduce the former 
we work in the massless theory but renormalize the associated mass operator by inserting it 
in a quark 2-point function at zero momentum insertion. This was the procedure used in the 
original three loop MS renormalization of [28, 29 . We then use the renormalization condition 
that there is no finite part at the subtraction point. In this computational setup we can still use 
the Mincer algorithm, |12l I13j . Thus we can deduce the renormalization constant and hence 
the three loop quark mass conversion function which is 



C m (a,a) 



l + C F 



-a 



4] a 



+ 



-18a 2 C j4 + 24a 2 C F 



288C(3)Cf + 57Cf + 332N f T F 



84aC A + 96aC F + 432C(3)Ca 
C F a 2 



1285C A 



24 



,,3^2 



71685a 2 C\ 



-\3l22a 6 C\ + 15552a 3 C A C F - 7776a 3 C F + 8748C(3)a 2 C A 

- 23328C(3)a 2 C A C F + 89424a 2 C A C F + 46656C (3) a 2 C F - 31104a 2 C£ 
+ 103032C (3) aCi - 357777aC A - 334368C(3)aC j4 C F + 573804aC A C F 

- 31104C(3)aC j4 iV/T F + l\34Q0aC A N f T F - 46656C(3)aC| - 30132aC| 

+ 62208( (3) aC F N f T F - 123120aC F N f T F + 3368844C(3)Cl - 466560C(5)C A 

- 6720046C A - 2493504C(3)C A C F + 155520C(5)C yi C , F + 2028348C A CV 

- 532224C(3)C , A A r /T F + 186624C(4)C j4 iV / T F + 3052384C A N f T F 

+ 451008C(3)C| + 933120C(5)C| - 2091096Cf, - 331776((3)C F N f T F 

- 186624C(4)C F N f T F + 958l76C F N f T F - 27648C(3)iV / 2 T| 
C F a 3 



240448 JV/Tj 



7776 



+ 0(a 4 



(3.10) 



Equipped with this and the result of |24^ 125] we find the minimal MOM quark mass anomalous 
dimension is 



7^ M0M («,«) 



3C F a + a 2 C A - 67C A - QC F + 8N f T F 



C F a 2 



+ 



- 3a 3 C\ + 24a 3 C A C F - 54C(3)a 2 C A + 411a 2 C A + 108a 2 C A C F 

- A8a 2 C A N f T F + 396C(3)aC A + \haC\ + 72aC A C F + 48aC A N f T F 
+ 5634C(3)Ci - 10095C A - 422A({3)C A C F + 244(7 A C F 

- 1152((3)C A N f T F + 3888C A N f T F - 3096C F + 1536( (3) C F N f T F 

C F a 3 



+ 736C F N f T F 



384iY/r| 



48 



+ 



1134a 4 C A CK(3) + 2835a 4 C A C F ({5) - 10125a 4 C A C F 

8586a 4 C A C| + QA8a 4 C 2 A C F N f T F - 2592a 4 C A C F 
35316C(3)a 3 C A C F - 1620C(5)a 3 C A C F + 45522a 3 C A C F 



- 15552C(3)a d C l A C F + 43740a 3 C A C F - 9072a 3 C%C F N f T F 
+ 31104C(3)a 3 C A C F - 7776a 3 C A C F - 485352C(3)a 2 C A C F 

- 28350C(5)a 2 C A C F + 893025a 2 C A C F + 167832( (3) a 2 C\C F 

+ 134784a 2 C A C F + 8h53QC{3)a 2 C 2 A C F N f T F - 307152a 2 ClC F N f T F 
+ 1 60704C (3) a 2 C A C F + 246888a 2 C A C F - 82944C(3)a 2 C A C F N f T F 

- 90720a 2 C A C F N f T F + 31 104a 2 C A C F NfT F - 4U72((3)a 2 C F N f T F 
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+ 799956C (3) aC\C F - 286740C(5)a(7 A C F + 295551aC A C F 

- 417744C(3)aCiC| + 223m8aC 2 A C F - 658368C(3)aC*|C F 7V/TF 
+ U54l6aClC F N f T F - 824256C(3)aC A Cf- + 432864aC A C 3 F 

+ 463104C(3)aCUCji7V / T F - 2505Q0aC A C F N f T F 

+ lQ5888((3)aC A C F NfT F - 46080aC A C F NJ?T F 

+ 248832C (3) aCpNfTp + 20736aC F N f T F - lW592((3)aC F NfT F 

+ 27M8aC 2 F NfT F + 16036470C(3)C A C F - 6334605C(5)C A C F 

- 10139319C A C F - 10029096C(3)ClC| + 3421440C(5)C A C| 

- 2188530<7 A C| - 15748128C(3)CiC F iV / r F + 2737152C(4)C A C F iV / T F 
+ 4147200C(5)Cl C F N f T F + 8403640C A C F N f T F + 2208384C(3)(7 A Cf. 
+ 6842880C(5)C* A Cf - 4669704<7 A Cf. + 7091712((3)C A CpN f T F 

- 2737l52((4)C A CpN f T F + l2MW0((5)C A CpN f T F 

- 22U048C A C F N f T F + 2405376({3)C A C F Nj?T F 

- 995328C(4)C A C F iV / 2 r| - 2128192C A C F N^T F - 1741824C(3)C> 

- 817128C£ + 4935168C(3)C*|7V / T F - 7464960C(5)Cf,AT / T F 

+ 3509568C F N f T F - 1327W4((3)C F NJ?T F + 995328C (4) C^iV/T 2 , 

- 605568C F Nj?T F + 147456C(3)C l FA r / 3 T| - 2048C F N?T F 

j(4) j(4) ,(4) 

+ 1244160C(3)^ - 165888-^ - 2488320C(3)-^iV / 
N F N F N F 

d {A) 
+ 331776-^ N f 

N F 



a 



i 



5184 



+ 0(a b ) (3.11) 



which involves the same rank 4 Casimirs as the /3-function. We have checked the three loop part 
by the direct evaluation of the anomalous dimension using the minimal MOM renormalization 
constant. Given that we are considering a relatively new scheme we have also renormalized 
the flavour non-singlet vector current. This is important since it is a conserved physical cur- 
rent and its anomalous dimension is zero to all orders in perturbation theory. This is true in 
all schemes but we have checked this explicitly to three loops by repeating the above quark 
mass operator renormalization but using the vector current, tp^tp, instead. With the minimal 
MOM quark wave function renormalization constants and isolating the Lorentz channel of the 
Green's function with the inserted current corresponding to the transverse part, we have checked 
that the vector current renormalization constant is unity to three loops in the minimal MOM 
scheme. Thus the Slavnov- Taylor identity has been checked to this loop order with the above 
renormalization. 

For more practical purposes it is useful to provide the explicit numerical expressions for 
SU(3). Thus we have 

C g (a, a) = 1 + [- 0.375000a 2 - 0.750000a + 0.5555567V/ - 7.041667 



+ 



+ 



0.210937a 4 - 0.562500a 3 - 0.625000a 2 N f + 0.950346a 2 - 0.4166670^/ 

- 7.437954a - 0.15432107V} 2 + 24.491 186/Y/ - 163.35991 lj a 2 

- 0.131836a 6 + 0.791016a 5 + 0.585937a 4 N f - 7.768301a 4 - 0.937500a 3 N f 

- 20.064612a 3 - 0.17361 la 2 Nf - 18.480594a 2 N f + 8.788748a 2 

- 2.423078a/V/ + 56.307562aiV/ - 530.662942a + 0.085734JV/ 
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C A (a,a) 



C c (a, a) 



C$(a,a) 



C m (a,a) 



1 + 
+ 



47.5818307V/ 2 + 1099.9356417V/ - 5176.895449J a 3 + 0(a 



0.750000a 2 + 1.500000a - 1.1111117V/ + 8.083333 



2.812500a 3 + 8.437500a 2 - 1.66666707// + 31.418274a + 1.2345687V/ 2 



53.9129287V/ + 256.103491] a 2 



+ 



19.411733a 4 + 81.359870a 3 



13.754707a 2 7V> + 272.349881a 2 



+ 6.929489aTV 2 - 254.788106a7V- + 1621.114903a - 1.3717427V, 



7 



+ 171.2676487V/ - 2601.1663737V- + 9357.562431 



7 



'/ 



a 3 + 0(a 4 



1 + 3.000000a + 
+ 



1.346529a 2 + 3.072567a - 5.9375007V f + 80.935770 



9.344565a 3 + 61.885373a 2 



7 



5.501305a7V f + 211.462123a + 8.7658777V 



7 



7 



431.8041367V/ +2945.691833] a 3 + 0(c 



1 - 1.333333aa + 



+ 



V 2 

.2 



-2.722222a 2 - 11.575317a + 2.3333337V/ - 25.464206 



16.906900a d - 72.363802a z + 23.739312a7V/ - 317.382214a 



- 6.4609057V/ 2 + 246.4426507V/ - 1489.980500 a 6 + O(a^) 
1 + [-1.333333a - 5.333333] a 



+ 
+ 



-1.222222a 2 - 6.888889a + 9.2222227V/ - 149.040228 
-11.953704a 3 - 44.682370a 2 + 14.024098a7V f - 269.395219a 



f 



11.7319307V/ 2 + 713.3336517V/ - 5598.952656 



a 3 + 0(a 4 



(3.12) 



Clearly it would appear that the series have large corrections at three loops. Though that for 
the quark wave function is best. 



4 SU(N C ). 

Although we have given the minimal MOM scheme results to as high a loop order as is possible 
for an arbitrary colour group, it is possible to provide the complete set at four loops for the case 
of SU(N C ). This is because the four loop MS anomalous dimensions of the gluon, ghost and 
quark are known for this colour group for an arbitrary linear covariant gauge fixing, [241 [26] . 
Using the electronically available data files associated with the latter article we have extended 
the various three loop minimal MOM results using the same method. This is also possible since 
we have the mapping for the gauge parameter between the two schemes at three loops. Thus 
for the gluon we have 



7^ M0M («,«) 



[3aTV c - 137V C + 47V/] 



6 



+ 



+ 



-6a 3 7V c 3 + 17a 2 7V c 3 - 8a 2 7V c 2 7V/ + 17aTV c 3 - 8aTV c 2 TV/ - 1707V C S 



+ 927V 2 N f - 247V, 



\f 



247V 



-165a 4 7V C 5 + 12a 4 7V c 4 7V/ + 54C(3)a 3 7V c 5 - 126a 3 7V C 5 - 72a 3 7V c 4 7V/ 

- 576C(3)a 2 7V c 5 + 1761a 2 7V c 5 + 72C(3)a 2 7V c 4 7V/ - 972a 2 7V c 4 7V/ 

+ 216a 2 7V c 2 7V/ - 774C(3)a7V c 5 + 102aTV c 5 - 144C(3)a7V c 4 7V/ - 588aTV c 4 7V/ 



13 



+ 



+ 288aN r 2 N f + 3456C(3)iV c 5 - 23032iV c 5 - 1080((3)N c 4 N f + 15500N c 4 N f 
- 1360A^ 3 A^ + 4224C(3)iV c 2 iV / - 47687V c 2 iV / - 768C(3)JV c iV f 2 + 73W c Nf 

1728C(3)a 5 iV c 7 - 7560C(5)a 5 iV c 7 - 23490a 5 iV c 7 + 1728a 5 N®N f 



+ 52389C(3)a 4 iV c 7 - 2250C(5)a 4 iV c 7 - 130590a 4 Nj - 1440C(3)a 4 iV c 6 iV j 



- 5040C(5)a 4 /V c 6 iV / + 1 0440a 4 N?N f + 3888C(3)a 4 iV c 5 + 2430C(5)a 4 iV c 5 

- 3888a 4 N 4 N f + 38214C(3)a 3 iV c 7 + 40410C(5)a 3 iV c 7 - 20727a 3 Nj 
+ 23616C(3)a 3 N?N f + U40((5)a 3 N^N f - 86328a 3 N?N f 

+ 1944C(3)a 3 /V c 5 + 53460C(5)a 3 /V c 5 + 3888a 3 N*N f - 937188C(3)a 2 iV c 7 
+ 656370C(5)a 2 iV c 7 + 1387395a 2 Nj + 324144C(3)a 2 iV c 6 iV / 

- 30240C(5)a 2 Af c 6 A( f - 1222596a 2 iV c 6 iV / - 18432C(3)a 2 iV c 5 iV / 2 

+ 93888a 2 N*Nf + 93312C(3)a 2 iV c 5 + 145800C(5)a 2 iV c 5 - 5832a 2 iV c 5 

- 319680C(3)a 2 iV c 4 Af / + 373104a 2 N 4 N f + 55296<C(3)a 2 iV c 3 iV / 2 

- 52992a 2 N*Nf + 5184a 2 N?N f - 3727014C(3)aiV c 7 

+ 1557630C(5)aiV c 7 + 1030995aiV c 7 + 1860768C(3)aiV c 6 iV / ' 

- 142560<(5)a/V c 6 /V / - 1 254696a N?N f - 400896C(3)aiV c 5 iV / 2 

+ 192288aiV c 5 iV / 2 + 872856<C(3)aiV c 5 - 461700C(5)aiV c 5 - 23328aiV c 5 
+ 24576C(3)aiV c 4 iV / 3 - 6144a7V c 4 iV f 3 - 383616C(3)aA^ 4 A(f 
+ 428544a N*N f + 82944C(3)a/V c 3 iV / 2 - 7948&aN?Nf + 7776aN 2 N f 
+ 8025711C(3)iV c 7 + 4451400C(5)iV c 7 - 27205691iV c 7 

- 5654160C(3)iV c 6 iV / - 4647600C(5)iY c 6 iV) + 23030340/V c 6 7V / 
+ 642816C(3)iV c 5 jV / 2 + 737280C(5)iV c 5 iV / 2 - 42469447V C 5 7V / 2 

- 2563488C(3)iV c 5 - 4949910C(5)iV c 5 + 142344iV c 5 - 24576C(3)iV c 4 iV / 3 
+ 1886727V c 4 iV / 3 + 8424000C(3)7V C 4 7V / - 2730240C(5)iV c 4 iV / 

- 745934AN?N f - 1810944C(3)iV c 3 iV / 2 + 552960C(5)iV c 3 iV / 2 
+ 2Z77728N* Nf + 110592C(3)iV c 2 iV / 3 - 165888AT c 2 J\r) 3 

- 1976832C(3)iV c 2 iV / + 3041280C(5)iV c 2 iV / - 326736/V c 2 iV / 

- 221184C(3)iV c iV / 2 - 552960C(5)iV c iV / 2 + 337536^^/ 

+ 59616JV ']207iW + O< ° 5 » C' 1 * 

where we have substituted the SU{N C ) values for Cf and Ca- Similarly, the ghost anomalous 
dimension is 



niMOM , 



, a 



7 C (a, a) = N c [a-3\- 



+ 



3a 2 N c - 3a N c - 3AN C + 47V) " 
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54C(3)a 3 iV c 3 - 45a 3 iV c 3 - 36C(3)a 2 /V c 3 + 216a 2 iV c 3 - 24a 2 N?N f 

+ 42C(3)aiV c 3 + 109aiV c 3 + 48C(3)aiV c 2 iV / - 76aN 2 N f + 564C(3)iV c 3 
- 51967V C 3 + 48C(3)iV c 2 iV) + 187QN 2 N f - 160N C N^ - 72N f 



a 

192 



"c 
3 



14 



+ 



3843C(3)a 4 iV c 5 + 2070C(5)a 4 iV c 5 - 8052a 4 iV c 5 + 72a 4 N c 4 N f 

- 3456C(3)a 4 iV c 3 + 1890C(5)a 4 iV c 3 + 8298C(3)a 3 iV c 5 - 6390C(5)a 3 iV c 5 

- 9411a 3 iV c 5 - 576((3)a 3 N c 4 N f + 192a 3 N c 4 N f + 10152C(3)a 3 iV c 3 

- 18900C(5)a 3 iV c 3 - 108C(3)a 2 iV c 5 - 32790C(5)a 2 iV c 5 + 73071a 2 N c 5 

- 2U8C(3)a 2 N c 4 N f + 3360C (5) a 2 N*N f - 32388a 2 N c 4 N f + 2880a 2 N 3 Nf 
+ 3888C(3)a 2 iV c 3 - 22680<C(5)a 2 iV c 3 + 648a 2 A^, 3 + 1296a 2 iV^iV/ 

- 40458C(3)aiV c 5 + 26790<C(5)aiV c 5 + 88801aiV c 5 + 1152C(3)aiV c 4 iV / 

- 12480C(5)aiV c 4 iV / - 42248aiV c 4 iV / + 3072C(3)aiV c 3 iV / 2 + 3136aN?N? 

- 1080C(3)aiV c 3 - 52380C(5)aiV c 3 + 2592aiV c 3 - 1728C(3)aiV c 2 7V/ 
+ 3600aA^ c 2 A^ + 310041C(3)iV c 5 + 192240C(5)iV c 5 - 1888893iV c 5 

- 91728( {3) N c 4 N f - 59040C(5)iV c 4 iV / + 9970687V C 4 7V/ + 13824C(3)iV c 3 iV / 2 

- 141984iV c 3 iV / 2 - 526176<C(3)iV c 3 + 549990C(5)iV c 3 + 14904iV c 3 

+ 3840iV c 2 iV / 3 + 54720C(3)iV c 2 iV / _ 34560((5)N 2 N f - W4928N 2 N f 



,4 



• 4608C(3)iV c iV / 2 + 18816iV c Af/ - 432JV,] ^^ + 0{a 5 ) . (4.2) 



Finally, the quark anomalous dimension is 

7 ^ M0M (a,a) = «[iV c 2 -l]^ 



+ 



+ 



+ 



3a 2 iV c 4 - 3a 2 iV c 2 + 6aiV c 4 - 6aN 2 + 22iV c 4 - 4N?N f - 19 N 2 

18C(3)a 3 iV c 6 + 15a 3 /V c 6 - 18C(3)a 3 iV c 4 - 3a 3 iV c 4 - 12a 3 N 2 

- 90C(3)a 2 iV c 6 + 105a 2 /V c 6 + 24a 2 iV c 5 iV / + 90C(3)a 2 iV c 4 

- 87a 2 N c 4 - 24a 2 N 3 N f - 18a 2 N 2 - 618C(3)aiV c 6 + 431aiV c 6 

+ 96C(3)aiV c 5 iV / - 32aN c 5 N f + 618C(3)aiV c 4 - 467aiV c 4 - 96((3)aN 3 N f 
+ 32aiY c 3 iY/ + 36aiV c 2 - 1182C(3)iV c 6 + 3231iV c 6 + 192((3)N?N f 

- 944iV c 5 JV/ + 32N*Nf + 894C(3)7V C 4 - 26377V C 4 - 192((3)N 3 N f 

i a 3 
+ 968N 3 N f - 32N 2 Nf + 288C(3)iV c 2 - 576/V c 2 - 24J%N f - 18j ——3 

330C(3)a 4 /V c 8 + 75C(5)a 4 iV c 8 - 151a 4 7V C 8 - 12a 4 NjN f - 420C(3)a 4 7V c 6 

- 45C(5)a 4 iV c 6 + 280a 4 7V C 6 + 12a A N^N f + 42C(3)a 4 /V c 4 - 30C(5)a 4 iV c 4 

- 105a 4 iV c 4 + 48C(3)a 4 /V c 2 - 24a 4 N 2 + 1116C(3)a 3 iV c 8 - 180C(5)a 3 iV c 

- 33a 3 /V c 8 - 192C(3)a 3 /V c 7 iV / + 184a 3 iV c 7 iV / - 1620C(3)a 3 iV c 6 
+ 180C(5)a 3 /V c 6 + 663a 3 iV c 6 + 192( (3) a 3 N c 5 N f - 184a 3 iV c 5 iV / 

+ 360C(3)a 3 /V c 4 - 558a 3 N 4 + 144C(3)a 3 iV c 2 - 72a 3 N 2 - 1662C(3)a 2 A^ 

- 192((3)N 3 N f + 3130C(5)a 2 /V c 8 - 5767a 2 /V c 8 - 384( (3) a 2 Nj N f 

- W0({5)a 2 NjN f + 2580a 2 Nj N f - 96a 2 N c 6 Nj i + 1518C(3)a 2 iV c 6 

- 4510C(5)a 2 7V c 6 + 6307a 2 iV c 6 + 576C(3)a 2 iV c 5 iV / + 160C(5)a 2 /V c 5 /V f 

- 2940a 2 iV c 5 iV / + 9§a 2 N 4 NJ + 144C(3)a 2 iV c 4 + 1380C(5)a 2 iV c 4 

- 558a 2 iV c 4 - 192((3)a 2 N 3 N f + 360a 2 N 3 N f + 18a 2 iV c 2 
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■s 



- 33680C(3)aTV c 8 + 17420C(5)aTV c 8 + U292aN c 8 + 7040((3)aN r 7 N f 

- 3200C (5) aNj TV/ - 3752aTV c 7 TV/ + 320aTV c 6 TV/ 2 + 49248C(3)aTV c 6 

- 37500C(5)aTV c 6 - 13642aTV c 6 - 9408C(3)aTV c 5 TV/ + 5760C(5)aTV c 5 TV/ 
+ 2632aTV c 5 TV/ - 192aTV c 4 TV/ 2 - 15568C(3)aTV c 4 + 20080C(5)aTV c 4 

- 122aTV c 4 + 2368C(3)aTV c 3 TV/ - 2560C(5)aTV c 3 TV/ + 1024aTV c 3 TV/ 

- 128a7V c 2 7V/ 2 - 528a7V c 2 + 96a7V c 7V/ - 172560C(3)TV C 8 + 61875C(5)TV C 
+ 252104TV C 8 + 43392C(3)TV C 7 TV/ - 12000C(5)TV C 7 TV/ - 118572TV C 7 7V/ 

V f 2 + 153927V C 6 7V/ 2 + 118002C(3)7V C 6 - 29085C(5)TV C 6 



1536C(3)A^ 6 A^ 2 + 15392TV C 6 TV/ 2 + 118002C(3)TV C 6 
225318TV C 6 - 640JV c 5 iV / 3 - 36864C(3)TV C 5 TV/ + 480C(5)TV C 5 TV/ 



+ 120868TV C 5 TV/ + 2304C(3)TV C 4 TV/ 2 - 16592TV C 4 TV/ 2 + 66462C(3)TV C 4 

- 55830C(5)TV C 4 - 19731TV C 4 + 640TV C 3 TV/ 3 - 6912C(3)TV C 3 TV/ 

+ 11520C(5)TV C 3 TV/ + 7152TV C 3 TV/ - 768C(3)TV C 2 TV/ 2 + 1200TV C 2 TV/ 2 

- 2304C(3)iV c 2 + 7680C(5)/V C 2 - 3974TV C 2 + 384C(3)TV C TV/ - 9448TV C TV/ 



9600C(3) + 15360C(5) - 3081] 



384TV 4 



+ 0(a 5 



(4.3) 



For practical purposes it is perhaps more appropriate to provide the explicit numerical values 
for all renormalization group functions at four loops for the SU(3) colour group. Thus 



/? mM0M (a,a) 



„,mMOM 



(a, a) 



[0.666667TV/ - 11.000000] 



+ 



+ 



- 2.250000a 3 - a 2 N f + 7.500000a 2 - aTV/ + 9.750000a 

+ 12.666667TV/ - 102.000000] a 3 

0.375000a 4 N f - 15.468750a 4 - 2.250000a 3 N* - 5.547924a 3 



V 
29. 170372a 2 N f + 151.166003a 2 



\f 



35.750000a7V/ + 174.652162a 



- 19.3833107V/ 2 + 625.3866707V) - 3040.482287 



+ 



2. 250000a 5 TV, 



/ 
61. 925145a 3 N 



114.265701a 5 + 4.816305a 4 TV/ 



298.616620a 4 



/ 



37.364446a 3 + 43.033474a 2 Nf 



+ 5540.175086a 2 + 3.385091aTV 



/ 



83.916446aTV/ 2 



1495.393156a 2 N f 
152.826933aTV/ 



- 1111.191853a + 27.4926407V) 3 - 1625.4022437V/ + 24423.3305507V/ 

- 100541.058601] a 5 + 0(a 6 ) 
[1.500000a + 0.6666677V/ - 6.500000] a 



+ 



- 2.250000a 3 - a 2 7V/ + 6.375000a 2 - aN f + 6.375000a 
+ 11.1666677V/ - 63.750000] a 2 



0.375000a 4 TV, 



/ 



15.468750a 4 



2. 250000a 3 TV 



26.920372a 2 TV/ + 100.182677a 2 
14.3833107V 2 + 444.8524147V. - 1769.783563 



</ - 5.727087a d 
22.784256a7V/ - 77.661754a 



7 



'/ 



a 



+ 



2.250000a 5 N f - 114.265701a 5 + 3.972555a 4 N f - 270.114333a 4 



'/ 



72.936261a 3 TV/ + 287.225188a 3 + 31. 783474 a z 7V, 



2AT f 



1126.940395a 2 TV/ + 3789.309068a 2 + 3.385091aTV/ 3 
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7, 



mMOM 



(a, a) 



mMOM 



(a, a) 



7, 



mMOM 



(a, a) 



124.724674aA/ + 1081.645997aiV/ - 6926.344667a + 22.492640A/ 



- 1141.450868iV/ + 14846.203053^/ - 54060.225189 
[0.750000a - 2.250000] a 

1.687500a + 0.750000iV/ - 19.125000 



0(a 5 



+ 
+ 



1.687500a 2 
2.799995a 3 

+ 22.427774a - 
0.421 875a 4 N f 



+ 5.625000a 2 N? 



1.125000a 2 iV/ + 24.289587a 2 - 0.857872aiV/ 



2.500000AT/ + 90.267128iV> - 635.349362 



V 



'/ 



26.892596a 4 - 2.931942a 3 N f - 121.006552a 3 

185.757176a 2 N f + 648.958841a 2 + 13.337341aJV, 



\f 



f 



- 314.266897aA?/ + 1090.833574a + 2.500000TV/ - 241.975687A/ 
+ 4788.563749^/ - 23240.416706] a 4 + 0(a 5 ) 



1.333333aa 

+ 



3.000000a 2 + 6.000000a - 1.333333A?/ + 22.333333 



9.492589a 3 + 2.000000a 2 iV/ - 0.296280a 2 + 6.949789a7V/ 
- 78.967792a + 0.888889ivf - 59.211534iV> + 459.481285 



7 



+ 



0.750000a 4 N f + 61.650284a 4 - 2.924683a 3 N f + 210.615979a 3 



2.00000a 2 A/ + 121. 134099a 2 N f 



869.566016a 2 + 6.962963aiV/ 



+ 77.834748aiV/ - 1553.524496a - 4.444444A/ + 281.560058A/ 



4934.050066iV/ + 20300.851595] a 4 + 0(a b 



4.0a + 



a 2 + 1.333333JVf - 69.666667 



+ 



+ 



1.916667a J 



/ 



2.000000a 2 iV/ + 98.522232a 2 + 2.000000aiV/ 



+ 130.753633a - 2.666667A/ + 152.122739iV/ - 1520.596003 



0.750000a 4 N f 

.2 at-2 



36.419738a 4 - 10.500000a 3 N f + 127.577470a 



+ 6.000000a z TV/ 
- 551.791827a7V/ 



345.848075a 2 N f + 3588.203465a 2 + 20.550022aiV f 2 



5040.515124a + 5.632797A/ 



'/ 



156.909331iV v 



f 



- 1073.781658iV/ - 9337.969739] a 4 + 0{a 5 ) (4.4) 

We have checked that the Landau gauge expression for the /3- function agrees with that of [6]. 

One interesting consequence of these expressions is that we can provide the anomalous di- 
mension of a particular dimension two operator which is 



O 



1 Aa Aaii 



(4.5) 



It is known, (30 |, I3 H [32] . that O has a novel renormalization property. In the Landau gauge the 
anomalous dimension of O is the sum of the gluon and ghost anomalous dimensions. Moreover, 
in an arbitrary linear and nonlinear covariant gauge there is a simple generalization of this 
Slavnov- Taylor identity which was established in |33| . This was based on the observation given 
in [33j which were explicit three loop MS computations. The operator is of interest as it was an 
attempt to have a gluon mass term in the Lagrangian which while not gauge invariant is in fact 
BRST invariant, [35] . It has seen renewed interest more recently, since it is believed to be the 
origin of dimension two power corrections in the running of an effective coupling constant in the 
low energy limit, [HI [9]. In [26] the four loop MS Landau gauge result was given. However, the 
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arbitrary a MS expression for a linear covariant gauge fixing was not recorded. As the gluon 
and ghost propagators are examined in the minimal MOM scheme, [36], it is worth providing 
the renormalization for the operators explicitly. Though for reasons of space we provide the 
SU(3) expressions 



MS 



Tff (a, a) 



SU(3) 



[27c 


n + 8N f - 105] — 




+ 


108a 2 + 567a + 548iV} - 4041 


a 2 
48 


+ 


12393a 3 + 4374C(3)a 2 + 52488a 2 



22356aiV / + 17496C(3)a 
+ 268272a - 12080iV/ - 28512C(3)7V / + 437304^ + 13122C(3) 

a 3 

- 20413891 

J 1728 

+ [- 3011499C(3)a 4 + 4133430C(5)a 4 + 8030664a 4 + 20824614C(3)a 3 

+ 708588C(4)a 3 - 10431990<C(5)a 3 + 39936807a 3 - 2939328C(3)a 2 iV / 
+ 314928C(4)a 2 iV / - 8669268a 2 iY/ + 93612348C(3)a 2 

- 3779136C(4)a 2 - 18305190C(5)a 2 + 159478227a 2 

+ 3359232C(3)aAT / 2 - 2796768aiV f 2 - 60046272C(3)a./V / 

- 11337408C (4) aN f - 1271881440^/ + 612180666<C(3)a 

- 43696260C(4)a - 513923130C(5)a + 1146415923a + 497664C(3)iV / 3 

- 846464JV/ - 19558656C(3)iV / 2 - 6158592C(4)iV / 2 - 107248896AT/ 

- 289945440C(3)iY/ + 104451120C(4)iV> + 313061760C(5)7V/ 

+ 20828935807V/ + 72636831C(3) - 46766808C(4) - 1908463680C(5) 

- 7232776173] ^^ + 0(a 5 ) (4.6) 
for non-zero a. Equipped with this then the equivalent minimal MOM scheme expression is 

^ MOM (a,a)\ su{3) = [27a + 8N f -W5]^ 



+ 
+ 



+ 



-108a 3 - 48a 2 iV/ + 387a 2 - 48aJV/ + 225a + 572iV) - 3978 

648a 4 N f - 26730a 4 - 3888a 3 N f + 21870C(3)a 3 - 31347a 3 

+ 3888C(3)a 2 iV> - 53136a 2 iV> - 102060C(3)a 2 + 337770a 2 

- 3888C(3)aiY/ - 36180aiY/ - 115182C(3)a + 43011a 

- 1536C(3)iV / 2 - 273287V/ - 29088C(3)iY/ + 959652]^ 

a 3 
+ 696924C(3) - 4993812] — — 

1728 

93312a 5 iV/ + 279936C(3)a 5 - 1224720C(5)a 5 - 3805380a 5 

- 77760C(3)a 4 iV> - 272160C(5)a 4 iV / + 557928a 4 N f 
+ 11078613C(3)a 4 + 1341360C(5)a 4 - 27025488a 4 

+ 1135296C(3)a 3 iV / + 77760C(5)a 3 iY/ - 4591728a 3 N f 
+ 13097214C(3)a 3 + 1319490C(5)a 3 - 10218393a 3 

- 221184C(3)a 2 iV / 2 + 1817280a 2 Nf + 14990832C(3)a 2 iV / 



48 



^The full expression for SU(N C ) is given in the attached data file. 
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- 816480C(5)a 2 iV / - 71613396a 2 N f - 149908644C(3)a 2 
+ 83215350C(5)a 2 + 277974261a 2 + 147456C(3)aiV/ 

- 36864a7V/ - 6801408C(3)aiV / 2 + 3556224aiV/ 

+ 98413056C(3)aiY/ - 10730880C (5) aN f - 75346200aiV) 

- 617646222C(3)a + 259312590C(5)a + 231547167a 

- 73728C(3)iV/ + 1125120AT/ + 8977920C(3)AT / 2 

+ 14254080C(5)iV / 2 - 82945888^/ - 276910992C(3)A f 

- 280604160C(5)iV / + 1438122060^/ + 1437422031C(3) 

a 4 
+ 816720570C(5) - 5780555523] — — + 0(a 5 ) (4.7) 

for SU(3). 

5 Discussion. 

We have provided all the renormalization group functions in QCD in the minimal momentum 
subtraction scheme introduced in [6]. To do this we have explicitly renormalized the theory and 
applied the renormalization prescription given in [6] to define the scheme. While [6] concentrated 
on the /3-function the other renormalization group functions are required for other problems 
such as the infrared structure of propagators and therefore we have provided that information. 
Currently the results are known at four loops for the SU(N C ) colour group and at three loops 
for a general group. One feature which differs from [6] rests in the renormalization of the gauge 
parameter. In [6j a was renormalized in the MS way whereas here we have chosen to follow 
a fuller approach and renormalize the gauge parameter according to the same criterion as all 
the 2-point functions. While this differs from [6] both sets of results agree in the Landau gauge 
which is the main gauge of interest for practical studies of the infrared dynamics of the gluon 
and ghost. 



Acknowledgements. We thank Dr A. Sternbeck for valuable discussions. 
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